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1. INTRODUCTION 


The Data Noising Program (IMP) Mill provide a neM program which simulates the 
effect of various noise sources on flight data. The program will read a 
Space Shuttle Functional Simulator (SSFS) format file, perform various noising 
operations on it, and rewrite a file in the SSFS format. The first noise 
function in^lenented was the time skew function, since this function deter* 
mines the first and last output data record. Following implonentation of the 
time ska< function, up to six additional data noising functions may be imple* 
mented; these are noise, location, bias, alignment, quantization, and sample 
rate. Three of tiiese possible six functions have been inq)1emented at the 
time of this documentation, and are described herein (time skew, noise, and 
quantization). 

The task to design, check out, and document the Data Noising Program was 
performed for the EX32 organization at the Johnson Space Center. 
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2. PR06RAN DESCRIPTION 


2.1 TIME SKEW FUNCTION 

2.1.1 PURPOSE 

The ONP introduces a time shift, ei Jier in a positive (right) direction or 
negative (left) direction, for a given measuraient, with respect to a normal 
(unskewed) measurement. 

2.1.2 INPUT DATA 

For a givei measurement, the skew time can be input positively or negatively. 
Hcmever, the magnitude of skew time can not be greater than the magnitude of 
absolute value of 5 times the different time between two records (at) positively 
or negatively. 

2.1.3 ANALYTICAL DESCRIPTION 
a. Positive skew time 
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Equation (1) is also applied for negative skew. The "lost" skew data will not 
be computed. Therefore, that data will not be written on the output file. 

2.2 NOISE FUNCTION 
2.2.1 PURPOSE 

This function is used to change a measurement by adding the noise. It includes 
the frequency, amplitude, and phase. 
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2.2.2 INPUT DATA 

For a given measurement* the following Information must be Input when the 
noise function Is applied 

a. Number of frequencies 

b. For each frequency: 

Frequency In Hertz: FREQ 
Amplitude of frequency: AMPL 
Phase angle In degree: PHASE 

c. Random number* which Is a 5-7 digit Integer where the last digit Is 
odd: NRANDM 

d. Type of random number: 

0 designates uniform 

1 designates normal: NTYPE 

e. Amplitude of standard deviation: STANDN 

2.2.3 ANALYTICAL DESCRIPTION 

Given a measurement (whether or not It has been applied by another function) for 
each number of frequencies. It will be applied by the following equation: 

Y = Y + AMPL(J) * SIN[FREQ(J) * TLOCAL + PHASE(J)] 

where 

Y = the measurement 

TLOCAL = local time, which Is zero from start of the output record. Time Is 
Incremented by At. 

In the case where the Input random number of this given measurement Is 0, 

Y will be the final value after noise function Is applied; otherwise* It will 
have the random noise added. 

Y^ = Y + STANDN * X 

where X = ZOR (NTYPE) 
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2.3 QUANTIZATION FUNCTION 

2.3.1 PURPOSE 

This function Is used to change a measurement from a continuous form to one 
with discrete steps to simulate digital data. 

2.3.2 INPUT DATA 

For a given measurement, the following Information must be Input wh«fi the 
quantization function Is applied. 

a. Number of bits of quantization: IBIT 

b. The upper and lower limits of range for which the measurement Is to be 
quantized: QMAX, QMIN 

2.3.3 ANALYTICAL DESCRIPTION 

For a given iiteasurement Y: If Y is less than QMIN, let Y = QMIN or If 

greater than QMAX, let Y = QMAX. Then Y Is the final value returned from 
the quantization function. Otherwise, two additional parameters must be 
defined as follows: 

a. DELTA - Number of steps Into which measurement range Is divided. 

DELTA * (QMAX - QMIN)/BITS 
where BITS * FLOAT (2**IBIT) 

b. INCR - Number of Increments of size DELTA Into which a measurement value 

will be divided. 

INCR = INT ((Y-(QMIN + DELTA/2.)) /DELTA ) 
and * QMIN + FLOAT(INCR + 1) * DELTA 

where 

Y = the measurement that has been applied by another function previously, or 

Y - the measurement from the Input data file. 

Y. - final value after it has been quantized. 

J 
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2.4 INPUT AND OUTPUT RECORDS FORMAT 


Input 

ftecord 


Time 
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N’Measureaients 


— P* 
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Output 

Record 


Time 
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K'Heasurements 
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K ■ (fH-K) Neasurements 


Each Input record consists of time and N-measurements. The output record will 
consist of time and M-measurements , where 

M = N + K 

where K is the number of measurements that are modified. 

(1 < K < 100) 
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3. PROGRAM USAGE 


The main program and Its subroutine subprograms reside In the secured file 
EX32-L72330*0APSPT. To execute the program, assuming that the user has 
assigned the Input data file as well as the output data file, the user 
enters: 

0XQT EX32-L72330*DAPSPT.RUN 

The following message Is displayed: 

PLEASE ENTER THE INPUT AND OUTPUT LOGICAL UNITS 

Key In tmo desired units (which must not be 5 and 6) that the user has 
assigned to the Input and output data files. The program then will read 
the header record. The last word of the header record Is the number of 
words In each subsequent record. Therefore, the following Is then displayed: 

THE NUMBER OF WORDS IN EACH DATA RECORD IS XXX 

The program will ask for the following Information on the display: 

PLEASE ENTER START AND END TIME IN MILLISECOND INTEGER 

The user will key In two legitimate times, and the following n^ssage will be 
displayed: 

PLEASE ENTER CHANNEL NUMBER 

The user enters any desirable channel number, and the program will then display 
the following function menu: 

PLEASE ENTER THE NUMERICAL VALUE REPRESENTING THE FOLLOWING FUNCTIONS 

1 - SKEW 2 - NOISE 

3 - LOCATION 4 - BIAS 

5 - ALIGN. 6 - QUANT. 

7 - SAMPLE RATE 

NOTE: Only skew, noise, and quantization are presently available. 
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B. If "1" Is entered, >the message 

PLEASE ENTER SKEW TIME IN MILLISECOND INTEGER 

appears on the screen, the program will always let the user know the current 
maximum or minimum value of skew time. The program then lists the function 
menu 

b. If "2“ Is entered, the following message Is displayed: 

PLEASE ENTER NUMBER OF FREQUENCIES 

Key In the legitimate value for number of frequencies, then the following 
message Is displayed 

PLEASE ENTER VALUES OF FREQUENCY IN HERTZ 

AMPLITUDE 
PHASE IN DEGREE 

to request the user to enter each set of frequency, amplitude and phase. 
Enter one at a time. The program will convert Internally Input values Into 
radians and print them out. After that, the message following Is displayed: 

PLEASE ENTER THE VALUES REPRESENTING THE VARIABLES AS 

* RANDOM NUMBER WHERE IT IS 5-7 DIGIT INTEGER, LAST DIGIT IS ODD 

* TYPE OF RANDOM WHERE 0 IS UNIFORM AND 1 IS NORMAL 

* AMPLITUDE OF STANDARD DEVIATION. 

Enter one set of those values; the program will print them out and return 
the user to the function menu: 

c. If "6“ Is entered, the message Is then displayed: 

PLEASE ENTER NUMBER OF BITS 

Key In the legitimate value, then the program will ask: 

PLEASE ENTER THE UPPER AND LOWER LIMITS OF RANGE FOR WHICH THIS CHANNEL IS 
TO BE QUANTIZED. 

Enter two values; the program then will print them out and Immediately 
return the user to the function menu. 
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In order to get out the function menu, the user simply enters "0", the 
program will ask for another channel number. However, If ”0" Is entered 
when the program requests the channel number, the program will then display 
the message: 

TOTAL NUMBERS OF CHANNEL ARE XXX 

From this point on, the output records are printed. Each output record 
consists of the number of channels In each Input record and the total 
number of modified channels. 
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4. DIAGRAM 
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Because the time skew function dictates the outcome of the output records, the 
above chart shows how the output file looks, depending upon the Input skew 
time. Five possible cases may occur: 

1 - No skew time 

2 - Input positive skew time Is less than maximum positive skew time 

3 - Maximum positive skew time 

4 - Maximum negative skew time 

5 - Magnitude of the absolute value of Input negative skew time Is less than 

magnitude of the absolute value of maximum negative. 
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